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Abstract 

Wall published a paper in 1960 on the Fibonacci sequence where 
he derived many results concerning the period and prime power di¬ 
visibility modulo m. His periodicity results have been generalized 
before to second order linear recurrences. Here we study the se¬ 
quences generated by such recurrences, with starting values of {0,1}: 
Among other things, we (a) derive new prime power divisibility re¬ 
sults not only for the general recurrence, but also interestingly for 
the squares of such sequences, (b) derive the period by new methods, 
(c) establish many novel identities modulo a prime power, (d) show 
derivations involving powers of matrices generated by these general 
recurrences, etc. 


1 Introduction 

Fibonacci sequences have been extensively studied in the past. While this 
prior work is highly interesting, we shall describe in this introduction only 
what is relevant to the main themes of this paper, which includes the exten¬ 
sion of the prime power divisibility and other results of Wall. After Wall [3], 
Robinson [4] derived periodicity using the Fibonacci matrix method (for a 
history of this method see [18]). General recurrences and over finite fields 
were studied by many including [5, 6]. Then came a breakthrough result of 
Matijasevich [8] that uses a Fibonacci identity (extended to general recur- 
rrences in Theorem 9). On a diffferent note, Knuth and Wilf [11] studied 
the interesting prime power divisibility of generalized binomial coefficients 
generalizing the well-known work of Kummer. The power of the logarithmic 
bounding implied by Wall’s prime power result was used in [12] to show that 
the Fibonacci sequences written in any base is not context-free (our result 
in the current paper generalizes this result as well in a strong sense). Gen¬ 
eral recurrences were also presented in a nice exposition in[19] and shown 
to share many of the well known basic properties of Fibonacci and Lucas 
sequences. Periodicity was studied in [16] in the context of a splitting field 
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characterization. Prime power divisibility modulo m was further studied 
in [13,14,15, 17] and many related results were proved, but none of these 
extend the prime power results of Wall to general recurrences. 

Here we study sequences generated by general second order linear recur¬ 
rences (with starting values of 0 and 1), defined by the equation e(n)=Ae(n- 
l)+Be(n-2) where A and B are integers. For our known Fibonacci sequence 
A and B are 1 with starting values 0 and 1. Among other things, we gener¬ 
alize the prime power divisibility results of Wall on Fibonacci recurrences 
modulo m, show new derivations involving powers of matrices generated by 
these recurrences, prove many novel identities, study the starting values, 
etc. 


2 Arithmetic Index Progressions in Recurrences 

Wall proved in 1960 [3, Theorem 3] that the indices of elements in the Fi¬ 
bonacci sequence that are zero modulo m form an arithmetic progression. 
We will extend Wall’s Method to obtain the same result for general recur¬ 
rences (this will be needed later). 

Theorem 1: For the general sequence e(n) = Ae(n — 1) -I- Be(n — 2) the 
terms e(n) = 0(mod m) forms a simple arithmetic progression except in 
the case for which gcd{B,m) ^ 1. 

Proof: We prove by induction that the following identity (well known for 
the Fibonacci sequence since B=1 for Fibonacci sequence) holds for the 
general recurrence e{n + t) = e(n -I- l)e(<) -I- Be{n)e{t — 1) (this like many 
other results is equally easy to prove with hbonacci matrices [19], but we 
reserve this for attacking some significant problems later). 

Obviously for e(0) = 0 and e(l) = 1 and e(2) = b the relation trivially holds 
true for f = 1 and t = 2. Say it holds for all t <= k. Now for t = A: -|- 1 we 
have that 

e(n -I- l)e{k + 1) + Be(n)e{k) 

= e(n -I- l)[Ae{k) -I- Be{k — 1)] -I- e{n)[B * Ae{k — 1) + B * Be{k — 2)] 

= A[e{n+\)e{k) + B*e{n)e{k — 1)] -|- B[e{n-\-l)e{k — 1) + B * e{n)e{k — 2)] 

= Ae{n -I- fc) -I- Be{n -I- /c — 1) 

= e(n -I- fc -I- 1) 

Hence the relation is true for e(n) using Wall’s ]3] arguments on setting 
e(i) = 0{mod m) and e(j) = 0(mod m) and finding that e{i + j) = 
0{mod m). But for the case where gcd{B,m) ^ 1, then the sequence goes 
into a cycle where e(n) does not equal 0 for all n. 
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This is because for the period to exist we need e(n)=0 and e(n+l)=l for 
some n. Now it means that B.e(n-l)=l(mod m) or e(n — 1) = _B“^(mod 
m) but B~^{mod m) does not exist.Hence the period does not exist. We 
can in fact state that 0 will also not appear, otherwise one can go backwards 
if gcd(a,m) = l and we get back 0 and 1 but this is also not possible since 
it goes into a cycle where 0 and 1 do not exist consecutively. But the fact 
that 0 will exist on the periodic cycle of e(n)(mod m) depends on the fact 
that one can go backwards in the cycle . So in this case the lemma does 
not hold. 


3 Prime Power Divisibility in General Sequences 

Here u(n) is taken to be our normal Fibonacci sequence and v(n) is taken 
to be the Lucas sequence with starting values 2 and 1. Wall [3] studied 
prime power repetition and divisibility in Fibonacci Sequences : 

Wall’s Theorem: If u{n) = 0{mod p®), but ^ 0{mod is the first 

occurrence for r such that u{r) = 0{mod p®) then u(pn) = 0{mod p®“''^) but 
^ 0{mod p®"*"^) is the first occcurence for u{x) = 0{mod p®"*"^). 
ProofiAccording to Binet’s formula u{n) = (r" — s")/'\/(5) where r and s 
are conjugate surds obtained by solving for the n’th Fibonacci number as 
a function of n and v{n) = r" + s". Also u{an) = (r“" — s“")/-\/(5) . So 
solving for r“" and s“" and using the binomial theorem and the property 
that v(n)=v(n-l)+v(n+l), we get that 
gcd{v(n),e{n)) = 1 or 2 and 
u{an) = 2^““ * {°'Cj)j=odd^^^~^''^'^u{nyv{nY~^ 

So u{an) = 2^~°‘u{n){Ku{nY + av{n)°‘~^) . 

Hence for u(an) to be divisible by p®’*'^ p must divide a given that p® divides 
u(n). 

We now extend Wall’s theorem above to the general recurrences. 

Theorem 2: In sequence e (E) defined before, if e(n) = 0{mod p^), but 
0{mod p^+^j is the first occurrence for r such that e(r) = Qimod p^) then 
e{pn) = 0{mod p^+^) but 0{mod is the first occcurence for x such 

that e(x) = 0{mod p^+^). 

Proof: For e(n) = Ae{n — I) + Be{n — 2), by Binet’s formula, we have 
e(n) = k{r'^ — s") where r + s = A and rs = —B (unless specified otherwise 
we will assume e(0) = 0 and e(l) = 1). We need to prove the existence 


3 



of a sequence v(n) such that v(n) = k'{r'^ + s”) for a constant k'. Taking 

up from here, we can define v{n) = pe(n — 1) + qe(n + 1). We know that 

v{n) = Av{n— 1) + Bv{n— 2). Also we have e{n + 1) = Ae{n) + Be{n— 1). 

Solving these 2 equations we have Bv(n)=e(n+l)(pB+q)-A.q.e(n). 

Assuming that gcd{A, B) = 1 , we reach the conclusion that B must divide 

q for v(n) to be an integer for all n. So let q = xB where x is an integer. 

Substituting, we have 

v{n) = {p + x)e{n + 1) — Axe{n) 

v{n) = k[{p + x)(r”+^ — s”+^) — Ax(r” — s”)] 

v{n) = fc[r”((p + x)r — Ax) + s”(Ax — (p + a;)s))] 

Now the coefficients of r" and s" must be equal. So we get 
(p + a;)(r + s) = 2Ax or (p + x)A = 2Ax ot p = x 
or p = q/B 

Hence the sequence exists v{n) exists if we choose p and q such that 
q = pB. So choosing p = 1, q = B we have 
v{n) = 2e(n + 1) — Ae{n). 

Now since gcd{A, B) = 1 we must have that gcd{e{n), e(n +1)) = 1 . Hence 
we must have that gcd{v{n), e(n)) = 1 or 2. 

We have e(an) = k(r^" - s“”) 

r” = (1/2) * {e{n)/k + v{n)/k') and s” = (1/2) * {—e{n)/k + v{n)/k') 

So, e(an) = (fc/(2)“)[(e(n)/fc + v(n)/fc')“ — {v{n)/k' — e(n)/k)°‘] 

Now, k' = k*[{p + x)r — Ax] . 

Since p = x = l,k' = k * (2r — A) or k' = k{r — s) or k' = 1. Since By 
Binet’s formula k = l/(r — s) /c'must be 1. 

Hence, e(an) = (fc/(2)“)[(e(n)/A: + u(n))“ — (u(n) — e(n)/A:)“] Also 1/kmust 
be of the form \/b because r and s are the roots of an equation whose coef¬ 
ficients are real. 

e{an) = 2^““ * {°-Cj)j=oddk^^~^^e(n)^v(n)°'~^ . 

So e{an) = 2^““e(n)(A'e(n)^ -I- au(n)““^) . 

This formula along with the fact that gcd{v(n),e(n)) = 1 or 2 along with 
the result of section on arithmetic progressions leads to the conclusion that 
the law of repetion of primes also holds for the general recurrence relation 
with the same argument as in the case of Fibonacci sequence. However we 
must remember that this proof is only for starting values of 0 and 1. 

While it is obvious then that there is always a member E(n) in the 
general sequence which has more trailing zeros (when written in base m) 
than any given k, we also get the following from the above theorem (please 
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see [12]): 

Corollary: The number of trailing zeros of any member E(n) in the general 
sequence written in any base m, is <= k.log{n) for some fixed k. 

This exponential gap can be a powerful tool in studying the recognition 
complexity of these and other sequences. While Automata Theory is not 
an object of study in this paper, we observe that this gap (even a non-linear 
one will suffice) is sufficient to prove that every such generalized sequence 
written in base m is non-context-free using the same technique as in [12]. 


4 Alternative Analysis using Fibonacci Matri¬ 


ces 


Ek= 


and the matrix A as 



E(] 


= 1 


This section tries to analyze an approach using the Fibonacci matrices cul¬ 
minating in a very interesting theorem when this approach is matched with 
that of the previous one which follows Wall’s work. Let us define the Fi¬ 
bonacci matrix 

^ e{k + 1) e(fc) ^ 

^ e{k) e(fc-l)^ 

So we have that En = A" * An where An = | ^ ^ 

\^e(0) 

So we have again Ena = * Aq 

or Ena = {En * (Ai)“^)“ * Ai 
or Ena = {EnT * 

Now taking determinant on both sides we observe that 
given that e(0)=0 and e(l)=l and e(-l) = l since e() is the fibonacci se¬ 
quence in this case. Hence, we have that 

e{na -I- l)e(na — 1) — e(na)'^ = (e(n + l)e(n — 1) — e(n)^)“ * (1)^~“ 
or e{na + l)e{na — 1) — e{na)^ = (e(n-I- l)e(n — 1))“ — e(n)^“ — a* e(n)^ * K 
where K contains the rest of the terms. Hence, we have that 
e{na + l)e{na — 1) — e{na)‘^ = (e(n -I- l)e(n — 1) — e(n)^)“ * (1)^~“ 
or e{na + l)e{na — 1) — e(na)^ = (e(n -|- l)e(n — 1))“ — e(n)^“ — a * e(n)^ * K 
So assuming that we have e(n) is the first number for which e(n) = 0 (mod 
p) but l{mod p^). Also from the previous section we know that e(np) is 
the first term for which e(.) is divisible by p^ 

On taking both sides congruent to and putting a=p, 
e{np+ l)e{np— 1) — e{npY = e{n + l)e{n — ly — e(n)^^ -bp* e(n)^ * K, for 
some K. From the previous section we know that if e{n) = 0{mod p) then p^ 
divides e(np). Using that fact we find that since p divides e(n) , p^ divides 
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p*e{n)‘^*K — e{n)‘^P. Since the gcd{e{n),e{n+l)) a.nd gcd{e{n),e{n—l) = 1 
, we have the following novel result which has to be satisfied if the theorem 
in the previous section has to hold true. 

Theorem 4: e{np + l)e(np — 1) = (e(n + l)e(n — l))P{mod p^) if e(n) = 
0{mod p) 

In that case, p^ divides the left hand side as well. But since the left hand 
side is the square of an integer hence p'^ must divide the left hand side.Hence 
p^ divides e(np) and e(np) is the first term which is divisible by p^ but not 
p^. So in this case the 2 approaches coincide. 

This idea can also be similarly extended to p® from p by using the same 
arguments and hence the extension is 

Theorem 5: e{np + l)e(np — 1) = (e(n + l)e(n — l)Y{mod if e(n) = 
0{mod p®) 


5 Recurrence Spaces 

We can define the concept of a recurrence space induced by {A, B}, as also 
had been previously done in [19]. 

Lemma: The sequence e(n) = A.e{n — 1) + B.e(n — 2) forms a space for 
which any 1 member forms the basis. 

Proof:- Say e(n) = A.e{n — 1) -b B.e{n — 2) and we define v{n) = R.e{p) + 
S.e{q) . 

Then the sequence v obtained as below also belongs to the space. v{n) = 
R[Ae{p — 1) -b Be{p — 2)] -b S[Ae{q — 1) -b Be{q — 2)] 

= v(n) = A[Re{p — 1) -b Se{q — 1)] -b B[Re{p — 2) -b Se{q — 2)] 

=> v{n) = Av(n — 1) -b Bv(n — 2) 

Corollary: This also holds for multiterm general recurrence relation i.e 
v{n) = aie(p) -b a^eiq) + aze{r) -b.... then v(n) belongs to the class V which 
follows the property e(n) = aie{n— 1) -b a 2 e{n — 2) -b a^e{n — ?>) + .... An in¬ 
teresting point to note here will be that the recurrence space defined above 
can also be extended to Fibonacci Matrices. We have defined v(n) to be 
v(n)=Re(p)+Se(q). 

So writing it out in the matrix form using the defintion of Fibonacci ma¬ 
trices given in the previous section, we have 
Vn = RAPRo + SA<iEo 
or Vo = RAP-^Eo + SAi-'^Eo 
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or Vo = BEo where B = RAP-" + S'A'?"” 
and we have Vn = A"BEo 
or Vn = A"Vo 

So this implies that recurrence space is also applicable to the Fibonacci 
matrices. 

This recurrence space concept will allow us to relate general sequences with 
different starting numbers to sequences with starting values as 0 and 1, as 
shown in few of the subsequent theorems 


6 Extension to the General Recurrence 


For the general recurrence , the technique used above to analyse the fi- 
bonacci matrices does not apply. The problem lies in the fact that e(-l) 
depends on the number m since B.e(-l)=l(mod m). So to extend this result 
to the general recurrence we have to take a slightly different path. 

Again we can define, 

'e(2) e(l)\ 


En = A" ^ * El where Ei = 


e(l) e(0) 


and A = I 

VI V 

Then we have that Ean = A"" ^ * Ei or Ean * A = A"" * Ei or Ean * A = 
{En*A*{Ei)<^-l)Y*Ei 

So taking determinant on both sides, we have 
El = —1 and A = —b 
thus yielding the following result: 

Lemma: e(na + l)e(na—1) —e(np)^ = (—&)““^*(e(n + l)e(n—1) —e(n)^)“* 

(-l)l-a 


Since we know from the previous sections that the theorem from Section 
4 holds for the general recurrence as well, it gives rise to a very interesting 
result below because only if it is satisfied the above equation will be 0(mod 
^e+i) q’jjjg inference can be drawn from the previous 2 sections and the 
above equation. 

Theorem 6: {—b)P-^ * (e(n + l)e(n — 1))^ = e{np + l)e{np — l){mod 
where e(n) = 0{mod p®). 
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7 Periodicity : Extension of Wall’s work to 
general recurrences 

Let us first review Wall’s work on periodicity once more: 

Wall’s theorem on periodicity: If k{p‘^) ^ k{p) (the well known Wall’s 
question) then k{p^) = p^~^k{p) and if t is the largest number such that 
kip*) = k{p) then kip^) = p^~*k{p) where e>t for the Fibonacci sequence 
Proof:-In earlier section we extended WallaAZs theorem for the Arithmetc 
progression of the zeroes in the fibonacci sequence modulo a particular 
prime power and itaAZs higher power to the general recurrence. 

Taking it up from there,Wall then proved that 

u{an + 1) = 2~^{K * u(n)^ + a * u{n) * v{n)°'~^ + u(n)“) 

We now know that u(nx) are the terms which are equal to 0(mod p®) 
and u(pnx)=0 are the terms which are equal to 0(mod So setting 

u(pnx+l) = l(mod p®’*’^) to obtain fc(p®“''^) = pnx the equation above gives 
{v{nx)/2) = l{mod p®^^). 

So v{nx) = 2(mod p®) and thus u(nx+l)=l (mod p®). So we have that 
fc(p®“''^) = pA:(p®) and Wall’s question whether fc(p^) = k{p) remains unan¬ 
swered 

So we work along the same lines as above and prove the following : 
Theorem 7: If fc(p^) A k{p) then fc(p®) = p^~^k{p) and if t is the largest 
number such that fc(p‘) = fc(p) then k{p^) = p^~*k{p) where e>t for the 
general sequence 

Proof: From earlier sections, we have obtained v{n) = r" -I- s" and we also 
have u(n) = fc(r" — s") where k=I/(r-s) and r and s are the 2 conjugate 
surds which helps to relate u(n) and v(n) as a function of n. THe exact 
values of r and s can be calculated by using generating functions or using 
matrix theory to find the eigen values of the transforming matrix. So again 
we have that, 
e(an -I- 1) = fc(r“"“''^ — 

e{an + 1) = (fc/(2)“)[r * {e{n)/k -I- v{n))'^ — s * {v{n) — e(n)/fc)“] 

Since it is exactly the same form as the one obtained by Wall, the same con¬ 
clusions also follow i.e we have e(an-l-I) = 2““(iL*u(n)^-|-a*u(n)*u(n)““^-|- 
v{n)‘^) after similar simplification. Again setting e(pnx+I)=l(mod p®“''^)) 
we have that iv{nx)/2) = 2{mod p®) adnd thus we have e(nx+I)=l(mod 
p®) from v(n)=2e(n+I)-Ae(n). So the same theorem applies for the general 
recurrence as well. However the question of whether fc(p^) = fc(p) exists for 



some p still remains unanswered. 


8 Periodicity using Fibonacci matrices 


Periodicity is easy to prove modulo a prime p since only adjacent pairs 
can exist in the sequence, but it is possible to show a much smaller period 
by using Fermat’s little Theorem as has been done by many authors. Pe¬ 
riods for Fibonacci matrices were proved by Robinson [4], and for general 
recurrences by [16] among others. In this section we are going to prove 
an interesting theorem on periodicity for the general recurrence using the 
Fibonacci matrices. 

Theorem 8: For a sequence e(n)=ae(n-l)+be(n-2) the square numbers 
e(0)^,e(l)^,... generated by squaring the members of the general recur¬ 
rence satisfies the property that fc(p®) = if k{p) and 

if t is the highest number such that k{p*) = k{p) then A:(p®) = * k{p) 

where e>t 

Proof: By definition of the Fibonacci matrix for the general recurrence, we 

"e(2) e(l)\ 


have En+in = A" * Ei and Epn+in = * Ei where Ei = 

^ A B' 


e(l) e(0) 


and A = 


1 0 


So by substituting like before, we have that 

Ep„+1 = iEr,)P * iEi)^-P 

Taking determinant on both sides and observing that determinant of Ei=-1 

and 1-p is an even number, we have that 

e{pn + 2)e{pn) — e{pn + 1)^ = (e(n -I- 2)e(n) — e(n -I- 1)^)^ 

So if e(n) is the first term such that e(n)=0(mod p®) then pn is the first 
term for which e(pn)=0(mod So if we take congruent to p®’*'^ on 

both sides and set e(n -I- 1) = 1(tooc? p®) so that A:(p®) = n we have ,by 
using binomial expansion that e{pn + 1)^ = e(n -|- l)^P{mod p®’*’^) and 
hence e(pn -|- 1)^ = 1(tooc? p®“''^) This concludes the proof of the theorem 
since it means that fc(p®+^) = pn = p®“^ * k{p^). Again the possibilty of 
k{p^) = k{p) remains open. 
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9 Starting Points of Period in Generalized 
Recurrence 

Going back to sectionl, we stated that 0 may not the starting point of the 
periodic cycle in e(n)=ae(n-l)+be(n-2) if gcd{b,p) 1 . So the interesting 
objective of this section is to investigate the nature of the starting point. 

Say that for the series e(n)=ae(n-l)+be(n-2) (mod p), the first 3 terms are: 

0,1,a. So, for l,a to be in the periodic cycle we must have , 
a=a+bx (mod p) 
or bx=0(mod p) 

or x=tp/gcd(p,b) (mod p) for some integer t ^ 0 since 0 cannot be in the 
cycle 

So the next question is the existence of e(n-l) if e(n)=x, e(n+l)=l and 
e(n+2)=a. For that we must have , 
tp/gcd{p, b) + be{n — 1) = l('mod p) 
or be(n-l) = l-(tp/gcd(p,b))(mod p) 

For the solution to exist, we must have that 
gcd(b,p) I (l-tp/gcd(p,b)) 

Therefore, 

tp/gcd(p,b)=l(mod gcd(p,b)) 

or t = {p/gcd{p,b))~^{mod gcd(p,b)) 

This is condition for t. Hence e(0) should be 
[{p/gcd{p, b))-^(mod gcdijp, b))*{p/gcd{p, b))]{mod p). 

We can also further establish the starting point of a periodic cycle for any 
arbitrary v(0) and v(l) by making use of our results for recurrence spaces 
which states that v(n) lies in the recurrence space of E for which e(n) is the 
basis. So we can define v(n)=Re(n-l)+Se(n) and calculate R and S from 
the information given about v(n) and we can then derive the starting point 
of v(n) as well since we know the starting points of e(n). 

10 Related Divisibility Results 

In section 4 we proved that e{an) = 2^~°'e{n){Ke{n)'^ + av{n)°'~^) where 
gcd{v(n),e{n)) = 1. The following result is an extension of a result used 
in the proof of Hilbert’s tenth problem [8,9]) to well-known general recur¬ 
rences. 

Theorem 9: e(n)^ divides e{nm) if and only if e(n) divides m 

Proof: From the result in section 4 we can easily observe that since gcd(v(n),e(n))=l. 
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e(n)^ divides e(an) if and only if e(n) divides a. Hence the result follows 


The following is a extension to general recurrences of th result due to 
Hoggatt et al. [16] 

Theorem 10: divides e(ne(n)^) 

Proof: The proof is by induction on k. The case k=l holds because of 
Theoreml. Now say that e(n)^ divides e(ne(n)^“^) is true. 

We can write e(ne(n)^) = 2^““e(ne(n)^“^)(itre(ne(n)^“^)^+e(n)u(n)®(”)“^) 
by substituting n = ne{ne{n)^~^) and a = e(n) in the result from sec¬ 
tion 4. Since e(n)^ divides e(ne(n)^“^) it follows that divides 

e(ne(n)^“^)(i4re(ne(n)^“^)^ -I- e(n)u(n)®^”^“^) and hence the theorem. 

The following well known result for the Fibonacci sequence also extends 
naturally to a general recurrence starting in {0,1} (also see [19] for a dif¬ 
ferent approach) 

Theorem 11: If the sequence E has starting values of 0,1, then E(a) divides 
E(b) if and only if a divides b. 

Proof:- From Theorem 1, we have that for a prime p, the sequence of in¬ 
dices for which E(a) = 0{mod p) form an arithmetic progression. Since 
E(0)=0 and a divides b then a and b must be members of the same arith¬ 
metic progession. Hence p also divides E(b). This is true for all primes 
that divide E(a). Hence E(a) must divide E(b). For the converse, that is a 
divides b, that E(a) divides E(b) can be proved by using Binet’s formula. 
E{a) = (p“ — q°')/k and likewise for E(b). So if E(a) divides E(b) then a 
must divide b. 


11 Conclusions 

Practical applications of the Fibonacci and Lucas sequences are well-known 
in diverse domains including many in computer science, number theory 
(including primality testing), etc.. In fact the generalized sequences are as 
fascinating if not more (as pointed out in [19]), and perhaps more potential 
impact. What we have done in this paper is to try to explore some of the 
harder areas like prime power divisibility, which could have implications for 
many open problems in these areas. 
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